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I.  Introduction 


The  possibility  of  using  hydrogen  masers  as  frequency  standards  on  board  spacecraft  re¬ 
quires  careful  consideration  of  long-term  reliability  issues.  Since  many  of  the  likely  failure  modes 
of  a  maser  involve  either  the  atomic  hydrogen  source  system  or  vacuum  pump  problems  due  to 
the  hydrogen  load,  the  maser's  reliability  will  be  enhanced  by  efficient  use  of  its  hydrogen  supply. 

To  operate  the  maser,  an  excess  of  hydrogen  atoms  in  the  F=  1,  A/  =  0  hyperfine  state  must 
be  continuously  fed  [1]  into  the  maser  bulb.  An  rf  discharge  dissociator  breaks  hydrogen  mole¬ 
cules  into  atoms,  and  the  state  selector  (typically  a  hexapole  or  quadrupole  magnet)  focuses  those 
atoms  having  F=  1,  M  =  0  at  the  maser’s  bulb  entrance  orifice  and  defocuses  those  having  F= 0, 

A/  =  0,  thus  creating  the  population  inversion  required  for  maser  operation.  Clearly,  the  first  re¬ 
quirement  for  efficient  hydrogen  use  is  that  a  large  fraction  of  the  hydrogen  molecules  flowing 
into  the  dissociator  exits  as  atoms.  But  additionally  there  is  a  subtler  requirement,  caused  by  the 
fact  that  the  focusing  properties  of  the  state-selecting  magnet  depend  on  the  atomic  velocities; 
since  faster  atoms  will  be  deflected  less  than  slower  ones  during  their  magnetic  field  traverse,  the 
“focal  length”  of  the  magnet  will  be  longer  for  those  faster  atoms.  Thus,  only  atoms  within  a  fairly 
narrow  range  of  velocities  will  be  focused  at  the  bulb  entrance  orifice  by  a  given  state  selector  de¬ 
sign.  If  there  is  a  mismatch  between  the  velocity  distribution  of  the  atoms  coming  out  of  the  dis¬ 
sociator  and  the  velocity-dependent  transmission  of  the  state  selector,  the  hydrogen  use  efficiency 
of  the  maser  could  be  seriously  impaired. 

Relatively  little  is  known  of  the  velocity  distribution  of  atoms  effusing  out  of  rf  discharge  dis- 
sociators.  In  many  cases  the  tacit  assumption  is  made  that  the  atoms  will  be  in  thermal  equilibri¬ 
um  with  the  dissociator  wall,  but  that  is  not  necessarily  the  case.  The  threshold  energy  for  molec¬ 
ular  dissociation  by  collision  with  electrons  in  the  discharge  plasma  is  about  8.5  eV  [2].  Since  the 
molecular  binding  energy  is  only  about  4.7  eV,  each  atom  carries  away  approximately  2  eV  of  ex¬ 
cess  kinetic  energy.  Depending  on  dissociator  geometry  and  gas  density,  the  hydrogen  atoms 
may,  or  may  not,  undergo  enough  bulk  and  wall  collisions  to  thermalize  fully.  Velocity  distribu¬ 
tions  of  atoms  exiting  dissociators  followed  by  cryogenically  cooled  thermal  accommodators  have 
been  measured  [3,4],  but  these  results  are  not  applicable  to  our  problem  since  such  accommoda¬ 
tors  would  not  be  used  in  space-qualified  masers. 

We  have  determined  the  velocity  distributions  of  hydrogen  atoms  effusing  out  of  an  rf  dis¬ 
charge  dissociator  having  a  geometry  and  operating  parameters  resembling  those  of  a  maser  dis¬ 
sociator.  We  have  also  performed  realistic  calculations  of  the  velocity-dependent  transmission  of 
hexapole-magnet  state  selectors,  investigated  designs  yielding  good  matches  between  atomic  veloc¬ 
ity  distribution  and  magnet  transmission,  and  explored  the  consequences  of  possible  mismatches. 
This  report  discusses  the  relevant  techniques  and  presents  some  results. 
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II.  Velocity  Distributions 


Atomic  hydrogen  velocity  distributions  have  been  determined  using  a  magnetic  deflection 
technique  described  in  detail  elsewhere  (5).  Fig.  1  shows  the  experimental  arrangement.  Hydro¬ 
gen  gas  is  fed  through  a  temperature-  controlled  Pd-Ag  leak  [6]  into  a  cylindrical  double-walled 
Pyrex  bulb,  15  cm  long  and  1.9  cm  in  internal  diameter.  Compressed-air  flow  between  the  walls 
provides  cooling,  and  rf  power  is  inductively  coupled  to  the  discharge  by  an  external  25-turn  coil. 
The  hydrogen  beam  exits  the  dissociator  through  a  0.1-cm-long,  0.025-cm-wide  slit,  is  collimated 
by  a  second  slit,  0.025  cm  wide,  set  at  d  =  63.7  cm  away  from  the  source  slit,  and  then  travels  be¬ 
tween  the  polepieces  of  an  L=  11.4-cm-long  electromagnet  configured  in  the  “two-wire”  geometry 
described  by  Rabi  et  al.  [7],  After  traversing  a  D  =  71.3-cm  drift  space,  the  beam  is  detected  by  a 
quadrupole  mass  analyzer.  Our  dissociator  operates  over  a  wide  range  of  rf  power  levels  and  hy¬ 
drogen  pressures,  as  shown  in  Fig.  2;  dissociation  fractions  of  up  to  80%  can  be  achieved. 

H - d - 1 - L—h - D - H 


HYDROGEN  “TWO-WIRE”  QUADRUPOLE  MASS 

DISSOCIATOR  MAGNET  FILTER  AND  DETECTOR 


Fig.  1.  Schematic  view  of  the  experimental  arrangement.  ITie  hydrogen  dissociator  can  be 
displaced  transversely. 

Molecular  hydrogen,  having  no  magnetic  dipole  moment,  will  travel  through  the  magnet 
without  deflection,  but  hydrogen  atoms  will  be  deflected  by  the  inhomogeneous  magnetic  field  in 
opposite  directions  depending  on  the  sign  of  their  magnetic  moments.  The  dissociator  is  at¬ 
tached  to  the  rest  of  the  apparatus  by  flexible  vacuum  bellows,  and  can  be  displaced  transversely 
by  micrometer  screws;  in  this  way,  the  angular  distribution  of  atoms  deflected  by  the  field  can  be 
measured.  We  have  shown  previously  [5]  that  if  /oW  is  the  distribution  of  detected  atoms  as  a 
function  of  source  slit  position  at  zero  magnetic  field  and  tfY)  the  speed  distribution  of  the  atoms 
leaving  the  source  slit,  then  the  distribution  of  detected  atoms  when  the  magnetic  field  is  turned 
on,f{X),  will  be  given  by 

m  =  X  \l  MX  -  WT 2 )  g(V)dv  ( i ) 

FM 
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where 


0.492  n 
a  m 


BLd  (1  +  -~) 
L  +  d 


(2) 


B  is  the  intensity  of  the  magnetic  field,  a  is  a  magnet  geometry  parameter  (2a  is  the  separation  of 
the  “equivalent  wires”),  m  is  the  atomic  mass,  and  n  the  effective  magnetic  moment  of  an  atom 
having  hyperfine  quantum  numbers  F,M,  given  by  the  Breit-Rabi  formula  |8],  For  hydrogen, 

x  +  M 

^  =  ±  (1  +  2Mx  +  ^F)l/2  (/> 


where  no  is  the  Bohr  magneton;  the  (  +  )  sign  corresponds  to  F  =  0,  A/  =  0,  and  F=  1,  Af  =  -1;  the  (-) 
sign,  to  F=  1,  M  =  0,  and  F=  1,  M-  1.  jc  is  proportional  to  the  ratio  of  magnetic-to-hyperfine 
energies,  x  =  B/( 507  G). 


Fig.  2.  Atomic  hydrogen  beam  fraction  vs  rf  discharge  power.  Total  dissociator  bulb  pres¬ 
sures  are  indicated. 

To  determine  the  atomic  speed  distribution,  f^X)  is  measured  first,  and  then  f{X)  is  calcu¬ 
lated  using  Eq.  (1)  and  model  velocity  distributions.  The  results  of  the  calculations  are  compared 
with  the  measured  f[X),  and  the  model  distributions  are  adjusted  for  best  fit  to  the  experimental 
deflection  data.  Fig.  3  shows  the  measured  undeflected  beam  profile  f{£X),  as  well  as  the  mea¬ 
sured  deflected  beam  profile  f(X)  at  a  field  of  995  G.  Our  technique  to  determine  velocity 
distributions  has  been  previously  validated,  and  the  apparatus  has  been  calibrated,  using  an  effu¬ 
sive  rubidium  beam  [5].  The  validation  experiment  showed  that  the  left  side  peak  in  Fig.  3  should 
be  the  one  used  to  fit  the  calculated  deflection  profile  f{X),  since  the  atoms  deflected  into  the  right 
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side  peak  travel  very  close  to  the  convex  polepiece  of  the  magnet,  causing  the  constant-force  ap¬ 
proximation  used  in  the  derivation  of  Eq.  (2)  to  break  down. 


x  (cm) 

Fig.  3.  Detected  atomic  hydrogen  flux  vs  dissociator  slit  position.  Full  line:  undefiected 
beam  (zero  field).  Circles:  beam  defl  cted  by  a  995-G  field.  Both  sets  of  data  have  been 
normalized  to  unit  height.  Dissociator  pressure:  0.058  Torr. 

The  analysis  of  preliminary  atomic  hydrogen  deflection  data  showed  [5J  that  the  atomic  ve¬ 
locity  distributions  are  significantly  narrower  than  Maxwellians;  since  then,  we  have  obtained  and 
analyzed  additional,  higher  quality  deflection  data,  and  the  results  support  that  conclusion.  The 
first  model  speed  distribution  used  to  generate  calculated  beam  deflection  profiles  was  a  beam- 
Maxwellian: 

g(V)  =  (2/l/0XW)5  exp[-  (K/K0)2]  (4) 

where  the  most  probable  velocity  Vq  was  treated  as  a  free  parameter,  since  the  temperature  of  the 
gas  within  the  dissociator  might  be  higher  than  the  wall  temperature.  An  example  of  “best  fit”  to 
the  deflection  data  achievable  using  a  Maxwellian  distribution  is  shown  in  Fig.  4.  Recalling  that 
very  fast  atoms  are  not  deflected  much,  and  very  slow  atoms  undergo  large  deflections,  Fig.  4 
shows  that  the  low  and  high  speed  wings  of  the  beam-Maxwellian  distribution  are  both  too  high. 
We  then  attempted  to  reproduce  the  deflection  data  using  a  Gaussian  as  a  model  speed  distribu¬ 
tion. 


gfY)  =  l/(^r  a)  cxp[-  (V-  W202) 


(5) 
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x  (cm) 


Fig.  4.  Atomic  hydrogen  (lux  at  detector  vs  dissociator  slit  position.  Circles:  as  in  Fig.  3. 

Full  line:  calculated  for  a  beam-Maxwellian  velocity  distribution,  at  7=  350  K.  Both  sets  of 
data  have  been  normalized  to  unit  height. 

where  the  most  probable  velocity  Vq  and  the  width  a  are  both  treated  as  free  parameters.  This 
model  reproduced  acceptably  the  preliminary  data  [5],  and  fitted  the  new  data  better  than  the 
Maxwellians,  but  still  failed  to  provide  an  acceptable  match,  due  to  the  smaller  error  bars  in  the 
present  measurements.  To  obtain  good  fits  to  the  data,  we  have  to  allow  the  speed  distribution  to 
be  asymmetrical.  This  was  accomplished  by  using  as  a  model  distribution  two  half-Gaussians 
joined  smoothly  at  their  peaks: 

gf} 0  =  1/(72*  o)  exp[-  ( V-V0f/2a\  ]forF<  V0  (6a) 

g(V)  =  1/(72 X  a)  exp[-  (V-  V0)}/2oj  JforC  >  V0  (6b) 

where  a  =  (a,  +  o2  )/2.  The  most  probable  velocity  Vq  and  the  two  partial  widths  a,  and  o2 
are  treated  as  free  parameters.  With  this  model  distribution,  acceptable  matches  to  the  deflection 
data  were  obtained  over  the  whole  range  of  dissociator  pressures  we  have  explored.  Fig.  5  illus¬ 
trates  the  match  obtained  between  calculated  and  measured  deflection  profiles  for  the  same  data 
shown  in  Fig.  4,  and  Fig.  6  shows  the  speed  distributions  used  to  calculate  the  deflection  profiles 
in  Figs.  4  and  5.  It  is  apparent  that  at  low  dissociator  pressures  the  atomic  speed  distributions 
are  nonthermal  and  significantly  narrower  than  Maxwellians. 
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Fig.  5.  Atomic  hydrogen  flux  at  detector  vs  dissociator  slit  position.  Circles:  as  in  Fig.  3. 
Full  line:  calculated  for  the  asymmetric  velocity  distribution  defined  in  the  text.  Both  sets 
of  data  have  been  normalized  to  unit  height. 


Fig.  6.  Velocity  distributions  used  to  calculate  the  deflected  beam  profiles  shown  in  Figs.  4 
and  5.  Full  line:  asymmetric  distribution,  as  defined  in  text.  Dashed  line:  bcam-Maxwcll- 
ian.  Both  curves  have  been  normalized  to  unit  area. 
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In  the  range  of  dissociator  pressures  we  have  explored  (25  to  350  mTorr),  neither  the  peak 
velocity  nor  the  partial  width  of  the  distribution  on  the  high-speed  side  change  significantly.  The 
dependence  of  the  low-speed  side  partial  width  on  pressure  is  shown  in  Fig.  7.  At  the  high  end  of 
our  pressure  range,  the  overall  width  of  the  distribution  is  quite  close  to  the  Maxwellian  width. 

The  peak  speed  of  the  distributions  corresponds  to  a  kinetic  energy  of  about  0.075  eV,  indicating 
that  the  hydrogen  atoms  do  lose  most  of  their  excess  energy  before  leaving  the  dissociator.  Since, 
on  the  other  hand,  the  atoms  do  not  thermalize  fully,  that  energy  loss  must  take  place  in  just  a  few 
collisions,  requiring  a  relatively  high  average  energy  loss  per  collision.  Impact  vibrational  excita¬ 
tion  of  hydrogen  molecules  within  the  dissociator  bulb,  with  an  energy  loss  of  0.546  eV/collision,  is 
a  likely  mechanism  of  atomic  slowing-down. 


Fig.  7.  Partial  widths  of  the  best-fitting  asymmetric  velocity  distributions  vs  total  dissociator 
pressure.  Circles:  low  speed  side.  Triangles:  high  speed  side. 
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III.  Hexapole  Magnets  as  Velocity  Filters 


A  hexapole  magnet  can  be  used  as  a  state  selector  in  a  hydrogen  maser  because  the  inhomo¬ 
geneous  hexapolar  field  will  exert  a  radial  force,  pointed  towards  its  axis  on  atoms  having  a  nega¬ 
tive  effective  magnetic  moment,  and  outwards  on  atoms  having  a  positive  effective  magnetic  mo¬ 
ment.  In  this  way,  for  a  well  chosen  combination  of  geometry  and  magnetic  field  strength,  atoms 
having  F=  1  and  M  =  0  can  be  focused  at  the  entrance  orifice  of  the  maser  bulb,  while  those  atoms 
having  F  =  0,  A/  =  ()  will  be  defocused.  The  focusing  conditions  are  velocity-dependent,  and  so  the 
state  selector  will  also  act  as  a  velocity  filter. 

The  force  acting  on  an  atom  immersed  in  an  inhomogeneous  magnetic  field  B  is  given  by 
F  =  grad(F)  |9|.  'Hie  effective  magnetic  moment,  given  for  hydrogen  by  Eq.  (3),  will  depend  on  the 
atomic  quantum  numbers  F,  M  and  on  the  magnetic  field  strength  B.  The  magnetic  field  strength 
within  a  hexapole  magnet  of  bore  radius  ro  and  poletip  field  B q  is  given  by  B  =  B^r/ro)2  [10],  The 
radial  equation  of  motion  for  an  atom  within  the  magnet  bore  is  then 

r  =  (n/m)  x  2(Bn/rl)r  (7) 

where,  for  the  M  =  ()  states  of  hydrogen,  p  =  ±x/(l  +x2)xl 2  ,  and  x  =  (F()/507  G)(r/r())2.  This  equa¬ 
tion  can  be  rewritten  in  a  form  more  suitable  for  numerical  solution  as 


ih  q1 

+  ft  (>J)1/2 


(O'  Q  =  0 


(8) 


where  Q  =  r/r(),  p;i  =  /i(j/(507  G),  and  w  =  (2po B[)/rtir^2ft2.  The  solutions  to  this  equation  oscillate 
about  the  magnet  axis  when  p<0,  and  are  outward  bound  when  p>0. 

Fig.  8  shows  the  maser  state  selector  geometry:  a  hexapole  magnet  of  length  L2  has  its  en¬ 
trance  plane  at  a  distance  L\  from  the  source  exit  plane;  the  entrance  plane  to  the  maser  bulb  is 
at  a  distance  L3  from  the  exit  plane  of  the  magnet.  The  axial  velocity  of  the  atom  is  constant.  The 
radial  velocity  is  constant  in  regions  a  and  c,  and  the  radial  acceleration  in  region  b  is  given  by 
Eq.  (7).  Atomic  trajectories  leading  from  the  source  into  the  maser  bulb  are  calculated  by  solving 
the  equations  of  motion  in  regions  a,  b,  and  c,  and  matching  radial  positions  and  speeds  at  the 
boundary  planes.  The  velocity-dependent  transmission  of  the  state  selector  can  be  calculated  by 
integrating  the  atomic  trajectories  leading  from  source  to  bulb  over  starting  coordinates  and 
angles.  Fig.  9  shows  the  transmission  curves  obtained  for  a  given  choice  of  parameters  for  atoms 
in  both  A/  =  0  states.  The  transmission  data  have  been  normalized  to  the  solid  angle  subtended  by 
the  bulb  entrance  orifice,  weighted  by  the  effusive  beam  angular  distribution  (i.e.,  it  is  measured 
in  units  of  the  total  atomic  flux  per  unit  speed  that  would  enter  the  maser  bulb  in  the  absence  of 
the  magnet).  For  high  speeds,  both  curves  converge  asymptotically  to  0.25,  which  is  the  statistical 
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weight  of  each  of  those  states.  The  broad  transmission  peak  for  F=  1  atoms  at  3800  m/s  contains 
those  atoms  which  undergo  one  radial  oscillation  in  their  traverse  of  the  magnetic  field;  narrower 
peaks  at  lower  speeds  contain  atoms  which  undergo  more  radial  oscillations. 


Fig.  8.  Schematic  view  of  the  state  selector  geometry  (see  description  in 
text). 


Fig.  9.  Transmission  of  a  hexapole  magnet  (in  units  of  beam  flux  at  the 
maser  bulb  entrance  orifice  in  the  absence  of  the  magnet)  vs  atomic  speed. 

Full  line:  F=  1,  Af  =  0  hydrogen  atoms  (selected  state).  Dashed  line:  F=0, 

M-0  hydrogen  atoms  (rejected  state). 

To  obtain  the  velocity  distribution  of  the  transmitted  atoms,  the  velocity-dependent  transmis¬ 
sion  of  the  magnet  has  to  be  folded  with  the  velocity  distribution  of  the  atoms  entering  it.  Let  us 
assume  that  the  atoms  leaving  the  dissociator  do  so  with  the  velocity  distribution  g^V),  illustrated 
by  the  full  line  in  Fig.  6  (i.e.,  the  distribution  yielding  a  good  fit  to  our  beam  deflection  measure¬ 
ments).  The  full  line  in  Fig.  10  shows  the  transmitted  distribution  when  the  magnet  has  been  opti¬ 
mized  forgo(T).  If,  on  the  other  hand,  the  magnet  had  been  optimized  for  a  beam-Maxwellian 
distribution  at  wall  temperature,  the  transmitted  distribution  would  be  the  one  shown  by  the 
dashed  line  in  Fig.  10.  Integration  of  the  transmitted  distributions  over  atomic  speed  yields  the 
atomic  fluxes  into  the  dissociator.  This  procedure  shows  that  assuming  the  atomic  speed  distri- 
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bution  to  be  thermal  introduces  a  mismatch  between  state  selector  design  and  atomic  speed  dis¬ 
tribution  which  causes  the  waste  of  65%  of  the  F—  1,  A/  =  0  atoms  leaving  the  dissociator. 


Fig.  10.  Velocity  distribution  of  the  hydrogen  atoms  focused  at  the  maser 
bulb  entrance  orifice.  Full  line:  state  selector  optimized  for  the  actual 
incident  beam  velocity  distribution.  Dashed  line:  state  selector  optimized 
for  a  Maxwellian  beam  at  the  dissociator  wall  temperature. 
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IV.  Conclusions 


Our  studies  of  velocity  distributions  of  hydrogen  atoms  effusing  out  of  rf  discharge  have 
shown  that  while  the  atoms  lose  most  of  their  excess  kinetic  energy  rapidly,  they  may  not  thermal- 
ize  fully  before  exiting.  At  low  dissociator  pressures,  the  atomic  hydrogen  beam  velocity  distribu¬ 
tion  is  significantly  narrower  than  the  beam-Maxwellian  distribution  characteristic  of  thermal 
equilibrium. 

Our  analysis  of  the  focusing  properties  of  hexapole  magnet  state  selectors  shows  that  the 
efficiency  of  hydrogen  use  by  the  maser  can  be  increased  significantly  by  optimizing  the  state 
selector  design  for  the  actual  atomic  hydrogen  velocity  distribution.  This  finding  is  of  particular 
importance  for  the  design  of  space-qualified  hydrogen  masers,  where  the  maser  reliability  is 
enhanced  if  hydrogen  consumption  is  reduced. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for  national  security 
projects,  specializing  in  advanced  military  space  systems.  Providing  research  support,  the 
corporation’s  Laboratory  Operations  conducts  experimental  and  theoretical  investigations  that 
focus  on  the  application  of  scientific  and  technical  advances  to  such  systems.  Vital  to  the  success 
of  these  investigations  is  the  technical  staffs  wide-ranging  expertise  and  its  ability  to  stay  current 
with  new  developments.  This  expertise  is  enhanced  by  a  research  program  aimed  at  dealing  with 
the  many  problems  associated  with  rapidly  evolving  space  systems.  Contributing  their  capabilities 
to  the  research  effort  are  these  individual  laboratories: 

Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat  transfer 
and  flight  dynamics;  chemical  and  electric  propulsion,  propellant  chemistry,  chemical 
dynamics,  environmental  chemistry,  trace  detection;  spacecraft  structural  mechanics, 
contamination,  thermal  and  structural  control;  high  temperature  thermomechanics,  gas 
kinetics  and  radiation;  cw  and  pulsed  chemical  and  excimer  laser  development, 
including  chemical  kinetics,  spectroscopy,  optical  resonators,  beam  control,  atmos¬ 
pheric  propagation,  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions,  atmospheric 
optics,  light  scattering,  state-specific  chemical  reactions  and  radiative  signatures  of 
missile  plumes,  sensor  out-of-field-of-view  rejection,  applied  laser  spectroscopy,  laser 
chemistry,  laser  optoelectronics,  solar  cell  physics,  battery  electrochemistry,  space 
vacuum  and  radiation  effects  on  materials,  lubrication  and  surface  phenomena, 
thermionic  emission,  photosensitive  materials  and  detectors,  atomic  frequency  stand¬ 
ards,  and  environmental  chemistry. 

Electronics  Research  Laboratory:  Microelectronics,  solid-state  device  physics, 
compound  semiconductors,  radiation  hardening;  electro-optics,  quantum  electronics, 
solid-state  lasers,  optical  propagation  and  communications;  microwave  semiconductor 
devices,  microwave/millimeter  wave  measurements,  diagnostics  and  radiometry,  micro¬ 
wave/millimeter  wave  thermionic  devices;  atomic  time  and  frequency  standards; 
antennas,  rf  systems,  electromagnetic  propagation  phenomena,  space  communication 
systems. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metals,  alloys, 
ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  nondestructive 
evaluation,  component  failure  analysis  and  reliability;  fracture  mechanics  and  stress 
corrosion;  analysis  and  evaluation  of  materials  at  cryogenic  and  elevated  temperatures 
as  well  as  in  space  and  enemy-induced  environments. 

Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray  physics, 
wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric  and  ionospheric 
physics,  density  and  composition  of  the  upper  atmosphere,  remote  sensing  using 
atmospheric  radiation;  solar  physics,  infrared  astronomy,  infrared  signature  analysis; 
effects  of  solar  activity,  magnetic  storms  and  nuclear  explosions  on  the  earth’s 
atmosphere,  ionosphere  and  magnetosphere;  effects  of  electromagnetic  and  particulate 
radiations  on  space  systems;  space  instrumentation. 


